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Figure  3.  Same  as  Figure  I .  except  for 
the  months  after  the  monsoon. 


A  Three-way  Comparison  of  Sea  Level  from 
the  Central  Equatorial  Pacific,  1985-89 


Introduction 

Temporal  and  spatial  features  of  central 
equatorial  Pacific  Ocean  sea-level  variation 
appear  similar,  in  measurements  from  two 
very  different  systems  (one  in  the  ocean  and 
one  carried  on  a  satellite),  and  in  results  from 
a  numerical  model  of  the  region.  In  particu¬ 
lar,  there  is  an  interannual  cycle:  during  El 
Nino,  Kelvin  waves  appear  at  the  equator, 
and  the  sea-surface  ridge  associated  with  the 
equatorial  current  system  shifts  southward: 
in  non-El  Nino  years,  instability  waves  ap¬ 
pear  at  6°N,  and  the  ridge  shifts  to  the  north. 
This  three-way  comparison  gives  support 
both  to  measurement  systems  and  to  the  nu¬ 
merical  model. 

In  Situ  and  Satellite  Measurements 
^nd  a  Numerical  Model - —  •  -  - 

Coincidentally,  from  early  19S5  until 
late  1989  there  were  two  quite  different  in¬ 
strument  systems  recording  sea  level  across 
the  central  equatorial  Pacific  Ocean.  One 
was  an  in  situ  array  of  inverted  echo  sounders 
(Chaplin  and  Watts,  1984),  island  subsur¬ 
face-pressure  gauges,  and  sea-level  gauges, 
maintained  near  160°W.  The  other  was  the 
Geosat  altimeter.  We  consider  the  records 
from  these  two  systems,  comparing  them 
with  one  another  and  with  the  output  of  a  3'/:- 
laver  numerical  model  of  the  global  ocean. 
During  the  five-year  duration  of  this  three- 
way  intercomparison,  there  occurred  a  warm 
El  Nino  event  in  1986-87  and  a  cold  La  Nina 
event  in  1988-89. 

The  in  situ  array  contained  three  invert¬ 
ed  echo  sounders  at  10°N,  8°N,  and  6°N,  and 
island  subsurface-pressure  or  sea-level 
gauges  at  6°N,  4°N,  2°N,  0.4°S,  4°S.  and 
9°S.  Since  the  array  has  -2°  spacing  only 
between  10CN  and  near  the  equator,  we 
consider  here  only  data  at  stations  in  this 
interval.  The  longitudes  of  these  stations 
ranged  from  157.5°Wto  162.5°W.  The  data 
were  processed  and  converted  to  sea  level 
(or  equivalently  dynamic  height),  as  de¬ 
scribed  in  Wimbush  el  al.  ( 1 990).  Data  gaps 
of  - 1  year  each  at  8°N  ( 1 986-87)  and  1 0°N 
(1988-89)  were  appropriately  patched 
(Donohue  et  al.,  1992). 

The  Geosat  time  series  were  computed 
using  a  combination  of  crossover  and  col- 
linear  methods  (Miller  and  Cheney,  1990). 
After  applying  the  standard  corrections  pro¬ 


vided  in  the  Geosat  geophysical  data  records 
(Cheney  et  al.,  1991a,b),  the  geodetic  mis¬ 
sion  data  were  processed  into  crossover  dif¬ 
ferences  and  the  exact  repeat  mission  data 
were  processed  into  collinear  differences,  to 
remove  all  time-invariant  signals.  Subse¬ 
quently,  a  linear  trend  was  removed  from  the 
differences  along  each  satellite  arc  over  the 
range  40°S-40°N  to  eliminate  time-varying 
orbit  error,  and  the  results  were  binned  in 
areas  8°  in  longitude  by  1°  in  latitude.  The 
sample  interval  in  a  typical  8°  x  1°  area  is 
approximately  2  days  for  both  types  of  data. 
After  performing  a  least-squares  adjustment 
to  the  crossover  data  and  shifting  the  col¬ 
linear  data  to  have  the  same  mean  as  the 
crossover  data  in  a  1-year  overlap,  the  com¬ 
bined  data  were  low-pass  filtered,  with  half- 
— power  periods  95  and  30  days,  to  produce 
data  sets  for  analysis  with  two  different  de¬ 
grees  of  smoothness. 

Means  of  both  the  in  situ  and  altimeter 
time  series  at  each  latitude  w-ere  set  to  match 
0-2000  db  dynamic  heights  calculated,  for 
these  latitudes,  from  mean  climatology 
(Levitus.  1982)  averaged  over  the  longitude 
range  of  the  array. 

The  3'/: -layer  global  ocean  model  is  a 
descendent  of  a  model  by  Hurlburt  and 
Thompson  (1980),  w  ith  extended  capability 
(Wallcraft,  1991).  The  horizontal  resolution 
is  0.35°  in  longitude  and  0.25°  in  latitude. 
The  model  was  spun  up  from  rest  to  statistical 
equilibrium  using  the  Hellerman  and 
Rosenstein  (1983)  monthly  w  ind  stress  cli¬ 
matology,  then  integrated  from  1981  to  1991 
using  monthly  averaged  European  Centre  for 
Medium-range  Weather  Forecasting  1000 
mb  w  inds  w  ith  the  1981-91  mean  replaced  by 
the  Hellerman-Rosenstein  annual  mean.  The 
model  is  thermodynamic  (prognostic  density 
equation)  to  allow  varying  stratification,  but 
the  only  thermal  forcing  is  a  relaxation  to  the 
annual  mean  density  climatology  (Levitus. 
1982).  Additional  results  from  this  simula¬ 
tion  can  be  found  in  Metzger  et  al.  (1992). 
For  our  comparison,  time  series  of  model  sea 
surface  height  were  sampled  at  0.1 -month 
intervals.  These  time  series  cover  0.5°S  to 
I0°N  at  0.5°  intervals  along  160°W. 

Intercomparison  of  Results 

Figure  1  shows  plots  of  all  three  time- 
series  sets  contoured  in  time-latitude  space. 
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(c)  Model 

Figure  I .  Sea  surface  elevation  for  the  time  period  1 985 -89.  for  (a)  in  situ  data  (60- 
day  low-pass  filtered),  (b)  Gensar  dn<a  (95-day  low-pass  filtered),  (c)  model 
output  (30-day  low-pass  filtered).  In  each  case,  elevation  is  contoured  as  a 
function  of  time  and  northern  latitude  (in  degrees)  Contour  interval  is  5  cm 
Positive,  zero,  and  negative  elevations  are  shown  as  solid,  doited,  and  dashed 
lines,  respectively. 


In  this  figure,  the'  in  situ  and  model  output 
time  series  have  been  low-pass  filtered  in 
visual  accord  with  the  smoothness  of  the 
Geosat  data  set.  The  in  situ  data,  Geosat  data, 
and  model  output  all  show  the  following 
features.  The  height  of  the  zonal  sea-surface 
ridge  (Wyrtki,  1974)  associated  with  the 
North  Equatorial  Countercurrent  varies  sea¬ 
sonally,  typically  being  maximum  in  Octo¬ 
ber  and  minimum  in  April.  From  the  ridge  to 
the  trough  (usually  near  4°N  and  10°N,  re¬ 
spectively),  the  height  difference  has  an  Oc¬ 
tober  maximum  of  about  25  cm.  From  year 
to  year,  the  position  of  the  ridge  varies:  it  is 
farthest  south  ( ~2°N)  during  the  warm  event 
(principally  in  1986),  and  farthest  north  ( 
~6°N)  during  the  cold  event  (in  1988). 
Also,  at  the  very  end  of  1 989,  after  the  Geosat 
data  end,  there  is  another  northward  move¬ 
ment  of  the  ridge  to  6°N  visible  in  both  the 
in  situ  data  and  the  model  output. 

Complex  demodulation  amplitudes 
(Koopmans,  1974)  were  computed  and  aver¬ 
aged  over  the  1985-89  duration  of  the  time' 
series.  In  Figure  2  these  average  amplitudes 
are  contoured  as  functions  of  period  and 
latitude.  Note,  in  all  three  cases,  the  two 
high-amplitude  peaks  at  about  6°N. 

The  in  situ  data  show  these  two  peaks  at 
periods  of  45  days  and  1  year.  The  45-day 
peak  is  probably  an  equatorial  current  system 
instability  (Philander  et  ai,  1985),  a  west¬ 
ward  extension  of  waves  originally  seen  in 
satellite  thermal  images  of  the  eastern  equa¬ 
torial  Pacific  (Legeckis,  1977).  The  annual 
oscillation  has  a  node  near  8°N  and  another 
peak  to  the  north  of  10°N,  but  with  the 
opposite  phase  (see  Figure  1 ).  There  is  yet 
another  slight  peak,  at  80  days  period,  cen¬ 
tered  on  the  equator  and  extending  out  to 
about  3°  latitude,  consistent  with  the  form  of 
an  equatorial  Kelvin  wave  (Fu  et  ai.  1991). 
Looking  at  these  amplitudes  year  by  year  (not 
show  n),  it  is  apparent  that  during  the  1 986-87 
warm  event,  this  Kelvin-wave  peak  is  espe¬ 
cially  strong  while  the  instability-wave  peak 
is  absent.  Also  the  instability-wave  peak 
period  varies  from  year  to  year,  being  short¬ 
est  during  the  1988-89  cold  event.  In  a  plot 
similar  to  Figure  la,  but  of  unfiltered  data  (not 
shown),  it  is  clear  that  these  waves  are  stron¬ 
gest  near  year's  end,  and  essentially  absent  in 
spring  when  ridge  height  is  a  minimum  and 
consequently  the  current  system  is  weakest. 

The  Geosat  data  show  the  same  general 
spectral  and  spatial  character:  the  annual 
peak  has  a  very  similar  latitudinal  structure, 
the  instability-wave  peak  is  again  near6°N, 


and  a  Kelvin-wave  peak  appears  at  80  days 
period.  There  is  an  additional  peak  at  about 
1 00  days  period  near  8°N.  A  similar  peak  at 
9°N  in  the  eastern  Pacific  has  been  previ¬ 
ously  observed  with  inverted  echo  sounders 
(Miller  et  ai,  1985).  Looking  at  the  Geosat 
amplitudes  year  by  year  (not  shown),  the 
peak  appears  at  8°N  during  the  1986-87 
warm  event,  the  time  of  the  in  situ  data  gap  at 
this  latitude.  Hence,  though  not  visible  in  our 
in  situ  data,  this  peak,  of  unknown  origin,  is 
probably  real  nevertheless. 

The  numerical  model  output  shows  the 
same  general  features  as  the  in  situ  and 


Geosat  data  sets,  except  the  equatorial 
Kelvin-wave  peak  and  8-9°N  peak  are  not 
visible.  The  instability-wave  peak  here  has 
a  longer  period  (70  days)  and  is  centered  at 
6.5°N.  Absence  of  the  Kelvin  peak  and 
shifting  of  the  instability  wave  toward  longer 
periods  may  be  results  of  the  monthly  aw,  - 
aging  of  the  w  ind-field  input. 

While  there  are  differences  in  the  de¬ 
tails,  it  is  encouraging  to  see  broad  agree¬ 
ment  in  several  fundamental  characteristics 
of  these  time  series  from  very  different 
sources.  Of  particular  interest  are  the  El 
Nino  cycle  variation  of  the  high-pressure- 
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(b)  Geosat 


Fig ure 2.  Complex  demodulation  ampli - 
tudes  of  sea  level  averaged  over  the 
time  period  1985-89  (except  for  a 
few  missing  weeks  at  the  beginning 
of  1985  and  the  end  of  1989  in  (a) 
and  (b)-see  Figure  I),  for  (a)  in 
situ  data,  ( b )  Geosat  data  (30-day 
low-p^ss  filtered),  (c)  model  out¬ 
put.  In  each  case,  amplitude  is 
contoured  as  a  function  of  period 
(in  days)  and  northern  latitude  (in 
degrees ).  The  lowest-level  contour 
shown  is2-5  cm  and  is  the  left-most 
contour  in  each  panel.  Contour 
inten’al  is  05  cm. 


ridge  latitude,  and  the  6°N  instability-wave 
peak.  In  satellite  infra-red  images,  this  insta¬ 
bility  is  clearly  visible  to  the  east  as  a  wave¬ 
like  disturbance  of  the  thermal  front  near  2°N 
(Legeckis,  1977;  Miller  et  al.,  1985). 

From  the  data  presented  here,  the  pres¬ 
sure  field  perturbation  associated  with  the 
instability  appears  greatest  near  6°N.  This  is 
farther  from  the  equator  than  (a)  the 
instability’s  frontal  expression,  (b)  pressure 
maxima  of  available  Iow-meridional-mode. 
free  equatorial  waves  (Fu  et  al.,  1991),  (c) 
regions  in  which  small-amplitude  perturba¬ 
tions  of  the  mean  flow  are  barotropically 
unstable  (Lukas,  1987;  Luther  and  Johnson, 
1990),  and  (d)  latitudes  of  high  eddy  kinetic 
energy  production  (Luther  and  Johnson, 
1990),  all  four  of  these  being  within  4°  of  the 
equator.  But  a  maximum  in  eddy  potential 
energy  production  in  the  thermocline  (90  m 
depth)  at  about  7°N  (Luther  and  Johnson, 
1990)  suggests  a  baroclinic  instability  in  the 
North  Equatorial  Countercurrent,  which  may 
be  responsible  for  the  peak  that  we  see.  This 
instability-wave  peak  is  not  visible  inMitchum 
and  Lukas's  (1987)  composite  Pacific  sea- 
level  spectrum,  because  there  were  no  mea¬ 
surements  in  the  range  4-7°N.  Nevertheless, 
it  has  been  noted  previously  in  Geosat  data 
(Perigaud,  1990).  We  believe  that  reproduc¬ 
tion  of  these  features  is  an  important  measure 
of  success  for  equatorial  models. 
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Introduction 

The  importance  of  long  equatorial  waves 
in  the  development  of  El  Nino  has  been 
w  idely  recognized  since  the  seminal  works  of 
Wynki  ( 1 975)  and  McCreary  ( 1 976).  Equa¬ 
torial-ocean  linear  models,  forced  by  ob¬ 
served  winds,  have  been  quite  successful  in 
reproducing  the  major  sea-level  signature  of 
El  Nino  (see,  e.g.,  Busalacchi  et  al„  1983). 
This  event  is  the  result  of  an  intrinsic  cou¬ 
pling  between  tropical  ocean  and  atmosphere 
(El  Nino-Southern  Oscillation,  or  ENSO), 
and  both  simple  and  complex  coupled  ocean- 
atmosphere  models  have  been  developed  to 
understand  the  phenomenon  (McCreary  and 
Anderson.  1991;  Neelin era/..  1992).  Analy¬ 
ses  of  model  solutions  have  led  to  several 
interesting  ENSO  theories,  such  as  the  de¬ 
layed-action  oscillator  (Suarez  and  Schopf. 
1988;  Battisti  and  Hirst.  1989).  in  which 
equatorial  Rossby  and  reflected  Kelvin  waves 
are  important.  Unstable  coupled  ocean-at¬ 
mosphere  modes  therefore  may  be  at  the 
origin  of  ENSO.  According  to  Wakata  and 
Sarachik  ( 1 99 1  ).these  modes  aresensitivetothe 
choice  of  Rayleigh  friction. 

Dissipation  in  the  form  of  Rayleigh  fric¬ 
tion  was  introduced  by  Gill  (1980)  in  an 
equatorial  long-wave  model.  For  a  multiple¬ 
mode  linear  model,  friction  is  usually  taken 
as  a  vertical-mode  dependent  parameter  (if 
possible  based  on  physical  assumptions)  to 
allow  separation  into  vertical  modes 
(McCreary,  1981).  Gente/u/.(  1983) discuss 


the  equivalence  of  horizontal  and  vertical 
damping  with  Rayleigh  friction  and  its  verti¬ 
cal-mode  dependence.  Through  a  model  best 
fit  of  the  observed  semiannual  zonal  current 
oscillation  in  the  equatorial  Indian  Ocean,  the 
same  authors  found  a  first-vertical-mode  de¬ 
cay  time  of  two  years.  Without  a  better 
estimate  of  damping,  many  modelers  have 
used  a  similar  value  for  the  Ray  leigh-friction 
coefficient. 

The  purpose  of  our  study  is  to  estimate 
the  Rayleigh-friction  coefficient  in  an  equa¬ 
torial  long-wave  model,  constraining  sea- 
level  results  to  best  fit  three  independent 
observed  sea-level  data  sets  (tide  gauges, 
moorings  and  GEOSAT).  Following  the 
results  of  Busalacchi  and  Cane  ( 1 985).  where 
the  addition  of  the  third  and  fourth  vertical 
modes  does  not  add  any  constructive  infor¬ 
mation  to  their  model/dala  intercomparison, 
our  study  is  done  using  two  vertical  modes, 
and  without  any  a  priori  hy  pothesis  on  the 
dependence  of  the  friction  coefficient  on  mode 
number. 

Data 

The  longest  sea-level  data  set  ( 1 975-89) 
used  in  this  study  is  deduced  from  daily 
means  of  tide-gauge  measurements  at  14 
islands,  situated  mostly  in  the  western 
Pacific  (Figure  1).  Between  instrumental 
uncertainties,  island  and  barometric  effects, 
the  error  in  sea-level  estimates  is  of  the  order 
of  a  few  centimeters. 


Daily  mean,  surface  dynamic-height 
fields,  relative  to  500db.  are  derived  from  the 
temperature  sensors  of  nine  TOGA-TAO 
moorings  and  three  current  meter  moorings 
(McPhaden  and  Hayes.  1 990).  using  Lev  itus 
(1982)  mean  TS  relations.  Most  of  the  dy¬ 
namic-height  time  series  cover  the  1986-91 
period,  except  at  0°-l40°W  and  0°-l  I0°\V. 
where  they  start  in  November  1983.  The 
significance  of  these  open-ocean  time  series 
is  altered  by  several  gaps  and  technical  con¬ 
straints  (such  as  the  use  of  a  reference  level, 
mean  TS  relations,  and  inadequate  vertical 
temperature  sampling),  which  probably  re¬ 
sult  in  error  ranges  at  least  equal  to  those  of 
the  tide-gauge  measurements. 

The  GEOSAT  altimetric  sea-level  data, 
resulting  from  the  geodetic  and  exact-repeat 
missions,  cover  the  April  1985  -  October 
1989  period  over  most  of  the  tropical  Pacific 
(Figure  1).  There  is  a  gap  in  the  data  from  30 
September  1986  to  8  November  1986  be¬ 
tween  the  two  different  missions,  and  there  i-- 
serious  data  degradation  after  May  1989. 
This  data  set  w  as  built  on  a  8°  longitude  x  1 " 
latitude  grid  and  interpolated  to  daily  values. 
Detailed  information  about  data  combination 
within  the  two  missions,  their  corrections, 
and  their  processing  can  be  found  in  Miller 
and  Cheney  ( 1990) and  Cheney  etal.(  1991a). 
The  improved  water-vapor  and  orbit  cor¬ 
rections  result  in  an  rms  difference  between 
moniiu  v  GEOS  AT  and  island  *ea  level  of  5  to 
4  cm  (Cheney  eta!..  1991b). 
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